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1 Introduction 



The AdS/CFT correspondence, or more broadly the notion of holography, has pro- 
vided us with a unique window into the dynamics of certain strongly coupled field 
theories [H |2l |3]. Viewed differently, the boundary field theory can be seen as a 
realization of a consistent quantum theory of gravity in an asymptotically Anti de 
Sitter (AdS) space. The finite temperature is easily implemented in AdS/CFT by 
introducing a black hole event horizon into the bulk. A black hole horizon plays 
the role of an infrared cutoff beyond which thermal fluctuations prohibit probing the 
vacuum. A particularly simple limit of the AdS/CFT correspondence is to focus on 
the low momentum/frequency regime where one expects quantum field theories to 
exhibit universal behavior dictated by hydrodynamics. More precisely, field theories 
at finite temperature and in the low momentum and frequency regime possess cor- 
relation functions with a universal structure fixed by hydrodynamics. Reproducing 
this structure from a weakly curved gravity background in asymptotically AdS space 
can be regarded as a test for AdS/CFT. This is so because the field theory side of the 
duality is strongly coupled even in the hydrodynamic limit. Through the AdS/CFT 
dictionary, operators in the boundary theory including fluctuations in the conserved 
quantities are mapped into various bulk gravity modes. Exploring near-equilibrium 
phenomena in the boundary field theory is holographically dual to studying pertur- 
bation theory of certain black hole spacetimes. 

In recent years, much effort has been devoted to reproducing linearized hydro- 
dynamics and diffusion physics in the context of AdS/CFT. For instance, various 
transport coefficients have been calculated in different boundary theories and ther- 
mal ensembles [1]. Some attention has also been paid to the hydrodynamic derivative 
expansion, such as second order hydrodynamics, still linear in fluctuations [5]. 

Non-linear hydrodynamics, on the other hand, is host to numerous interesting 
phenomena associated with the non-linearities in fluid dynamics, like turbulence. In 
the context of holography, it would be extremely interesting to see if there exists a 
gravitational and/or stringy dual to, for instance, a turbulent flow. 

Reproducing non-linear terms in the hydrodynamics expansion from gravity has 
been discussed in [B]. Having said that it is important to emphasize that in the 
context of AdS/CFT, phenomena associated with non-linear hydrodynamics have 
remained largely unexplored. There are many interesting phenomena associated with 
non-linearities in hydrodynamics, among them the so-called long-time tails which are 
the focus of this contribution. 

Imagine = {p,j"') to be a conserved current, where a,b,... denote the spa- 
tial coordinates. Consider the zero-momentum self auto-correlator of namely 
f d'^~^x{j"'(t,x)j^{0,0)), where d is the field theory spacetime dimension. Recall the 
constitutive relations for in the hydrodynamic regime 



where D is the diffusion constant. The first term detects spatial inhomogeneities and 
is diffusive, while the second term is convective, proportional to the local flow velocity 
u"". Omitted terms are suppressed at long- wavelengths by additional derivatives or are 
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cubic or higher in fluctuations. At zero-momentum, only the second non-linear term 
in dl]) contributes. It involves long-lived hydrodynamic fluctuations (diffusive modes 
and sound waves), which decay on diffusive time scales t ~ k~'^ at wavenumber k. This 
generates a power law tail in the correlator with an exponent which is a characteristic 
of diffusive physics 



The numerical prefactor in ([2]) is exactly calculable within the hydrodynamics frame- 
work and the topic will be reviewed below [7j. Notice that this phenomenon originates 
from statistical fluctuations around the equilibrium state so it should be suppressed 
by the entropy density s. For instance, in a deconflned gauge theory plasma s ~ N^, 
where Nc is the number of colors, the long-time tail will be an 0{1/N^) effect. Ac- 
cording to the AdS/CFT dictionary, this is mappeclfl to a one-loop effect in the bulk 
AdS quantum gravity. This was noted sometime ago by Kovtun and Yaffe |8j. The 
tail ([2]) cannot be seen in the classical gravity approximation. 

Conflrming the existence and reproducing this long-time tail from a one-loop bulk 
gravity computation is the subject of the present contribution. Interestingly, classical 
statistical fluctuations in the boundary theory will be mapped to quantum gravity 
fluctuations in the bulk. 

Quantum gravity is not a renormalizable theory so one might worry about the 
prospects of performing a one-loop computation within its framework. However, 
because we are dealing with an infrared phenomenon, our flnal result will not depend 
on an ultraviolet completion. The relevant one-loop integrals will be ultraviolet safe. 
For precisely the same reason, stringy corrections will play no role; pure Einstein 
gravity and Yang-Mills theory in Anti de Sitter space will suffice. 

It turns out that it is physically motivated to organize the one-loop computation 
in a causal manner such that it does not rely on the physics beyond the black hole 
event horizon. We come back to this point later in the text. 

The organization of this paper is as follows. We begin by reviewing the hydrody- 
namics predictions in section [2j We take time to lay the formal grounds of our work 
throughout the section [31 At the same time, we clarify why on an intuitive level our 
results are totally expected. The reader familiar with the background material could 
jump directly to the section 13. 7[ perhaps after reading subsection 13.11 Pre-requisites 
to compute the one-loop amplitude are discussed in section HI It includes material 
ranging from computation of the vertices to gauge-flxing. Section [5] gathers these 
ingredients. Our paper closes with concluding remarks in section [HI 

2 Review of hydrodynamic predictions 

We now sketch the calculation of the long-time tails in hydrodynamics. For details, 
we refer the reader to [7] which contains references to experiments (see also a nice and 

^Wherever a gravitational dual is available. 
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more recent account of the topic by Kovtun and Yaffe, in the context of relativistic 
fluids |H], which we will be following closely here). 

In a volume V, there are sV independent effective degrees of freedom where s 
is the entropy density. By the central-limit theorem, fluctuations of infrared modes 
are small and Gaussian-distributed. These fluctuations are fluctuations in conserved 
quantities. In particular, their equal-time correlators are set by thermodynamics, 
while their dynamics is governed by linearized hydrodynamics. At arbitrary time 
separation the full two-point functions, for conserved charges and local fluid velocity, 
are given by [TJ [Sj^ 

(3) 



V-'{piT,k)p{0,-k)) 
V-\u'{t, k)u^{0,-k)) 



T 

e+p 



k'y . 



+ 



k'y 



is the speed of sound. 



cos(fccsr) 



v/ + P) ^iid 



where S is the charge susceptibility, Cg 
7, = C/(e + p) + 2(f/-2)7,/(rf-l). 

Hydrodynamics at finite temperature is to be viewed as an effective statistical 
field theory. In > (2+1) space-time dimensions, hydrodynamics' canonical coupling 
constant has negative mass dimension. Therefore, the theory becomes weakly-coupled 
in the infrared [71 [9] . For this reason, a one-loop treatment becomes exact in the long- 
time limit Ifl Such a one- loop treatment amounts to approximating by its quadratic 
term 

' d'-'x {f{T, x)j\0)) ^ / d'-'x {p{T, xK{t, x) p{0)u\0)). (4) 



(5) 



This can be evaluated using the Gaussian two-point functions ([3]) 

{piT,x)u''iT,x)piO)u\0)) = (p(r,a;)p(0)) {u^ir, x)u\0)) . 
One finds [3 [8] 

d'-'x {r{T,x)j\0)) = (6a) 

e + p [47r(Z}+7^)|r|] — 

which holds up to terms suppressed by higher inverse powers of t. Here we have 
straightforwardly generalized the results of to space-time dimensions. 

Similar long-time tails are present in the stress tensor correlators. Consider the 
correlator J d'^~^x{1fy{t,x)t^y{Q)). There exists a second order non-linear term in the 
constitutive relations for the stress tensor analogous to ([T]) with t^"^ ^ [e+p)5u^5uy . 
One finds [H] 



d'^-\ (t"J'(r,x)t"^(0)) 



T2 



(d-l)(d+l) 



(47r7s|r|) 



1 

— + 



2d 



7r7r?kl)' 



(6b) 



^ In quantum field theory, the appropriate correlator is the symmetric correlator 
fc), p(0, — fc)}), which is manifestly real [S]. 

^ In (2+l)-dimensions, the canonical coupling becomes dimensionless and quantities run logarith- 
mically with scale. For instance, the effective viscosity becomes large logarithmically in the infrared. 
A Kubo formula for the viscosity does not exist. Nevertheless, the (3 function for the effective loop 
coupling turns out to be positive (at least for incompressible fluids), causing it to flow to a trivial 
fixed point in the infrared. This allows precise predictions to be made on long-time behavior [9]. 
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3 Long-Time tails from one-loop AdS gravity 



Our goal will be to reproduce (l6a|) and (l6bl) from a bulk one-loop gravity computation 
in AdS. We find the bulk diagrams which are dual to the hydrodynamic processes. 
We begin by describing the setup. 



3.1 Setup and notations 

We use lower case Latin letters i, j, . . . = 1 . . . d — 1 to denote the spatial coordinates. 
The Greek letters are reserved for both temporal and spatial coordinates and run over 
0,1 ... d — 1. Upper case Latin letters are used for the entire spacetime indices. The 
radial index will be referred to as "r". The superscript "(0)" will mean "background 
value" . The bulk action we consider has the following form 

^ f ^d+i^ r^fj. d{d-l) 1 



+ 7r^ / d'^Xy/^K. 



The background Ac/S'd+i-Schwarzschild black hole solves the equations of motion as- 
sociated to ([7]) 

where 

p{r) = coshd(27rTr), g(r) = tanh {2nTr). (9) 

In these coordinates, the black hole horizon is located at r = and T is its Hawking 
temperature. For simplicity, we will work in units where AnTi/d = 1. The chemical 
potential is turned off so the gauge fields vanish in the background|§ In particular for 
d = 4, this background is dual to strongly coupled Af = 4 super Yang-Mills theory at 
temperature T and vanishing R-charge chemical potential 

For the bulk gravitational perturbations Qf^i, — g^^} = h^y{r,x^), we adopt the 
Gaussian normal gauge (spacelike axial gauge) suitable in holographic settings 

hr^ = 0, hrr = 0. (10) 

A suitable gauge choice for the gauge fields is 

Ar = 0. (11) 
There will be still some residual gauge freedom rigid along the r direction^ 



^Even pure gravity in Schwarzschild Anti de Sitter space has enough structure to exhibit the 
long-time tail phenomenon. 

^ This residual gauge freedom will eventually make a come back and haunt us, however this will 
not be a major concern. 
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Figure 1: Schwinger-Keldysh time contour for the finite temperature path integraL 
3.2 Schwinger-Keldysh FormaUsm 

Our aim is to discuss real-time physics in the boundary CFT so we work directly in 
Minkowski signature. For this we need the so called Schwinger-Keldysh, or real-time 
formalism. This allows for a description of expectation values with prescribed initial 
states or density matrices which is appropriate for a thermal ensemble. A compre- 
hensive review of the method can be found in pi]] . In order to set our notations, here 
we include a lightning review of the subject mostly based on [TT]. 

To accommodate operators with real time arguments, the path integral time con- 
tour has to make an excursion to the real domain. Since only the initial state is 
prescribed, a doubled set of fields is required in the path integral: one is to imple- 
ment the forward time evolution up to the locations where operators are inserted and 
the second one to implement backward time evolution back to the time at which the 
initial state is prescribed. Averaging over a thermal ensemble of initial states is done 
by adding a vertical Euclidean part (with an extent /3 = 1/T) to the time contour, 
with periodic boundary conditions (antiperiodic for fermions). This leads to a time 
contour depicted in Fig. [TJ Denote the field values and sources on the two horizontal 
parts of the contour as 

(j)i(t,x) = (f)(t,x), (j)2{t,x) = (pit — ia), Ji(t,x) = J(t,x), J2(t,x) = J(t — id), 

(12) 

where the length of the vertical part of the contour is denoted by cr. Using this 
notation, one can write down a generating functional for the real-time correlators as 
follows 

Z[JuJ2]= f[V4>]e-x^{iS + i r dt f d'^xJi{x)(j)i{x) - I [ ^ dt [ d'^xJ2{x)(f)2{x)), 
J Jti J Jti J 

(13) 

where the Minkowski signature action S = J dtd'^xC[(j){x)] is evaluated along the 
contour. The fields Ji and J2 are background sources for 0i and its partner 02 
respectively. Due to the field doubling, the Schwinger-Keldysh propagator is a 2 x 2 
matrix given by second variation of the generating functional as usual 

^ / N 5^ In 2' , , , 

Gab{x-y) = a, 6 = 1,2. 14 

dJa{x)dJb{y) 

Note that a is arbitrary. We will choose a=0. At equilibrium, all these Green's 
functions are related to the retarded and advanced Green's function through Kubo- 
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Martin-Schwinger (KMS) relations [12] E 

Gil = -^Re Gr + (1 + 2nB)Im Gr, G22 = «Re Gr + (1 + 2?2B)Im Gr, 

G12 = 2nBlm Gr, G21 = 2(1 + nB)Im Gr, (15) 
where ub is the Planck factor 

and 

Gn{x-y) ^ tMx), </.(y)]) ^(a;°-y°). (16) 

3.3 Schwinger-Keldysh Formalism in AdS/CFT 

Thermal QFT in the Schwinger-Keldysh formalism and the Kruskal extension of a 
black hole spacetime are connected [13] . In a maximally extended eternal AdS black 
hole spacetime as depicted in Fig. |21 there are two causally disconnected regions R 
and L with two AdS asymptotic regions. They provide natural candidate regions 
for the doubled fields of the Schwinger-Keldysh formalism to dwell on. Standard 
AdS/CFT prescription [3] for computing boundary theory Euchdean amplitudes can 
be naturally extended to the case where the bulk geometry is maximally extended. 
In this case the bulk partition function generates the Schwinger-Keldysh correlation 
functions ( IT^ at finite temperature T|j This was accomplished by Herzog and Son 

mE 

There is a useful formulation of this, which was proposed both by Skenderis and 
van Rees [T7] and by Son and Teaney [18], in which only the R-quadrant needs to 
be considered explicitly. This is the formulation we shall use. Let us first recall that 
the AdS/CFT correspondence equates the full quantum gravity (string/M-theory) 
bulk path integral on AdS with prescribed Dirichlet boundary conditions with the 
generating functional of the boundary correlators 

^bdy[0bdy] = y"[l'0]0|i,dy=</,i,dyexp(z5'bulk[0]) = ^bulk[0bdy]- (17) 

Herzog and Son's approach identifies the path integral on the Kruskal plane with 
boundary conditions on the fields and their Keldysh partners in the quadrants R and 
L [11]. The boundary values of the bulk fields in quadrant R and L, ^J'^f, are identified 
as sources in the Schwinger-Keldysh formalism in the boundary theory. What it is 
important here to note is general choices of a can be accomodated, not only /3/2 
as was considered originally by Herzog and Son. The geometry which we shall use 

^Also known as fluctuation-disspation relations. 

^Provided the initial state for the Lorentzian black hole in the bulk is given by the Hartle-Hawking 
state [H]. 

method for calculating retarded functions had previously been proposed by Son and Starinets 
[15] in planar coordinates. 
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Figure 2: The maximally extended geometry of an AdS-Scliwartzshild black hole. 
The shaded region depicts the causal diamond associated to a causal boundary cor- 
relator. 

is precisely Herzog and Son's geometry, continued to cr = 0. This precisely brings 
quadrants L and R on top of each other. This continuation was described precisely 
by Son and Teaney [18j. 

Skenderis and van Rees' obtain the same geometry by considering directly the 
standard Schwinger-Keldysh contour depicted in Fig. [1] (see Fig. 2 of |1^ for a more 
precise description), convoluted onto the black hole exterior region. Independent of 
how one derives it, this geometry concretely realizes the physical picture that from 
the viewpoint of an exterior observer, a black hole is described by a density matrix 
rather than an S-matrix. This will allow us to study ordinary thermal field theory in 
the i?-quadrant. 

In the boundary theory, KMS conditions relate all correlation functions to the 
retarded propagators. The bulk theory in our limit is a "field theory" of interacting 
gravitons and gauge particles propagating in the AdS black hole spacetime. 

The KMS relations also hold non-perturbatively in the bulk black hole spacetime, 
for quantum fields in a fixed background. This has been confirmed long time ago by 
Gibbons and Perry [19] 

Gi2(r,r';p) = -mB(p)(GR(r,r';p)-G'^(r,r';p)), etc., (18) 

where Gab{f,r']p) is the partial Fourier transform of the bulk propagator function 
and r, r' are two arbitrary radii0 Relations f[T5]) will play a key role. 

Since here the metric itself fluctuates, we require an extension of the Gibbons and 
Perry's result. We note that to all orders, the background AdS black hole geometry 
admits a time-like Killing vector relative to which a "temperature" can be defined. 
The black hole Killing horizon will have constant surface gravity everywhere on the 
horizon. Although, bulk loop effects can renormalize background black hole's bare 

Note that we use the invariant combination — ^(Gr— Ga) in lieu of 2IniGR as in ([15]). The latter 
would be ambiguous in a mixed representation {r,uj), as it would depend on coordinate choice. 

^^In quantum field theory, the KMS relations hold non-perturbatively, however in quantum gravity 
it is not clear how to define "bulk operators" non-perturbatively. Therefore, the meaning of ([T8|) 
may be limited to a semiclassical expansion. 
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parameters (for instance, one expects Sth/th ~ ^{^/^T))-, they will not interfere 
with the KMS conditions. This is so because we keep the (Euclidean) length of the 
time-like Killing vector at the boundary fixed. This means that the computations 
are performed in the canonical ensemblejll Therefore, it is reasonable to expect the 
relation f lTSj) is true to all loop orders in quantum gravity. 

The retarded and advanced propagators are computable in the quadrant R. In 
practice, the KMS relations remove any concern about the interior. Put differently, 
to all orders, any treatment involving the full Kruskal geometry must agree with the 
predictions of an ordinary thermal field theory applied to the exterior region, provided 
it is such that the equilibrium condition is respected. 



3.4 Dirichlet conditions for quantum fields in AdS 

In AdS/CFT we need to impose Dirichlet boundary condition ( 1T7|) on the bulk path 
integral. We follow Symanzik's treatment of quantum field theories on spaces with 
boundaries. The boundary conditions are implemented by adding a boundary action 
|2U] . This method has been previously applied to certain interesting problems. The 
list includes studying Casimir effect [20] and critical phenomena in the presence of 
boundaries [21j. It will be manifest from this procedure that the boundary correlators 
defined by f|T7|l are in fact equivalent to the correlators of certain bulk operators, i.e., 
the momentum conjugates of the corresponding bulk modes0 

We illustrate the idea for a bulk scalar field with the following action 

5buik = -\ I rf''+'a;y^(V„0V> + 2V(0)), (19) 

where Tc denotes the position of the regularized boundary. Classically, the Dirichlet 
boundary condition 0(rc, x) = </>bdy(3^) is obtained from varying the action if a suitable 
boundary term is added to the action 

5bdy = / rf-^xv^^^^^i^ (0My(a;) - 0(r„a;)), (20) 

J r=rc V ^ 

where II,^ = —\/—gg'^^dr(l) is the conjugate momentum field for 0, 7^1, is the time-like 
induced metric on the boundary and 7 is its determinant. This suggests to use the 
path-integral 

^buik[</'bdy] = j [P</)]e'^f""['^l (21) 

with Sfuii = 5'buik+5'bdy as a replacement for (fTTl) . Now (pir^ x) is treated as dynamical 
inside the path integral. To justify ( 12T1) at the full quantum level, one must further 
check, since the boundary condition is obtained only as an equation of motion, that 
0(rc,a;) has no fluctuations around its mean (phdy^x). This can be readily verified 
at the free theory level -the combination (0(r, x)— ^bdy(2;)) converges to zero as a 
distribution when r — )■ Tc in any two-point functioi^fl -More generally, this is proved 



^^By equilibrium, the periodicity of the Euclidean time circle is the same everywhere. 
^^A similar observation was recently exploited by Iqbal and Liu |22j for one-point functions in a 
background field, in a context closely related to ours. See also [23] 
^^That is, when integrated against x-dependent test functions. 
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to all orders in perturbation theory around Eq. (5.14c) of [20] . 

Since the path integral depends on 0bdy only through S'bdy, from ( 1T7|) . (pi]) and 
( 12 ip . we see that boundary theory correlators (obtained as variations with respect to 
0bdy) are equal to the bulk theory correlators of II's 



Now the Dirichlet boundary conditions are implemented automatically in the path 
integral. It should be clear that the Schwinger-Keldysh indices ik are the same on 
both sides. This identity will be verified in a concrete example for two-point functions 
in section (14. 3p . The bulk loop expansion of the right-hand side of ( 122|) gives rise to 
so called Witten diagrams [5]. 

In our set up, we will be dealing with bulk tensor fields {Am, Iimn) rather than 
scalars. The construction here trivially generalizes to these cases. Since the cor- 
responding Yang-Mills and Einstein bulk actions are also quadratic in derivatives. 
In gravity, the conserved boundary currents and stress tensor (densities) are defined 

|25] as variations of the on-shell action due to variation of the boundary values of 
Qf,^ and 



where a boundary limit r — )■ is understood. In our gauge K^^, = ^l^^drQau ■ 
These operators, by definition, are the (radial) momentum conjugate to the metric 
and gauge fields^ They play the role of II^ in (|22|) . Note that the bulk equations of 
motion imply their conservation. 

The momentum conjugates II^ have (5-function contact singularities proportional 
to 6{ri—rj)6'^{xi—Xj) in their correlators. This forbids the r — )■ Tc limit to be taken 
directly inside (!22l) . In general, additional contact terms proportional to 6'^{xi—Xj) 
and its derivatives may be necessary to be added to the right-hand side of (122|) for the 
equality to hold. Within bulk perturbation theory, these contact term ambiguities 
as well as the bulk renormalization of (1221) are discussed by Symanzik [2Q|- These 
are distinct from the so called holographic renormalization issues, which arise when 
the limit Tc — ?■ oo is taken in keeping with the AdS/CFT prescription. Additional 
counter-terms, intrinsic to the boundary metric and matter, may be necessary when 
the regulator is taken to infinity [26]. This amounts to subtractions in (I23ap and 
(I23bp . It is worth emphasizing that these counter-terms will not affect our long-time 
tail discussion. 

^^Working with the tensor densities turns out to be more convenient than working with the tensors. 



{0'^{xi)---o'"MUy 



lim (m(ri,xi)---m"( 



'"n; ■^n j / bulk- 



(22) 




(23b) 



(23a) 
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3.5 Flow equations 



It will be useful to think of the operators in ( l23l) as extending into the bulk and not 
just as boundary operators. That is, we consider more generally 

t'u=t'uir), 3^=3^{r). (24) 

Note that at finite r, they lose their meaning as local current densities0 We derive 
simple radial flow equations for these currents in a very similar way to [22] . 

We only quote the flow equations we need. As expected, it turns out that the 
operators f l23ap and (I23b|) have simple radial profiles in the hydrodynamic limit. The 



following equation is easy to prove from Einstein's equations 

^gdr{^9K^) = ^'^R^^ + j,6'^^. (25) 

Here '^'^^R^^ is the Ricci tensor of the induced metric on the constant-r slices. Us- 
ing Yang-Mills' equations of motion and ( l25l) . one obtains the following radial flow 
equations for the traceless part of the stress tensor and the current fluctuations 

dr{5t^.) = ^5(v/^Wi?^g = + 0{u\ q') fiy^u, (26a) 
ovrG 

1 



driSf) = ^5{d,^F^^) = + 0{u', q'), (26b) 



9j+i 

where we have used the fact that '^'^'>R^^ oc evaluated on the background. Also in 
momentum space, 5^'^^R^J\ = 0{u!'^, g^), since fx and u are tangential to the constant-r 
hypersurface. It is worth noting that errors on the right-hand sides are not uniform 
in r. As we will see later, due to the near horizon oscillations in the field profiles 
oc r^*'^/^'^-^, the right-hand sides in fl26l) are small only up to exponentially small 
radius r ~ g-2vrT/<^^ 

In sum, in the hydrodynamic regime the radial profiles of the currents are simply 
constant over the bulk geometry expect for a an exponentially small region near the 
horizon. Clearly, it makes sense to organize the bulk computation in terms of these 
variables. 



3.6 What to compute 

Computing a one-loop fluctuation function in a black hole background directly would 
be very difficult: in following the past history of a fluctuation within perturbation 
theory, we would have to inevitably cross the black hole horizon in a finite proper 
time, rendering such a calculation daunting. However, we have argued above that the 
details of the interior region must be irrelevant at the end, as long as this region is 
treated in such a way that the KMS relations are respected. What we will compute 

^''There is no "local" stress tensor associated to the gravitational field. 
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(a) (b) (c) 

Figure 3: Feynman diagrams contributing to (a) current correlator (b) stress tensor 
correlator. Four-point vertices (c) will be found to have negligible effects. Wavy lines 
are bulk Yang-Mills fields and double lines are gravitons. 

is the zero-momentum retarded two-point functions of the bulk operators and t^y, 
at one-loop 

WG^^(r) = Jd'^-^x ^^^f{T,x,r^ oo) f{0,r^ oo))r, 

= / d'^-'x (t\(r, x,r^oo) t%(0, r ^ oo))r . (27) 

It should be completely obvious from f l22]) that the boundary thermal field theory 
inherits the KMS relations from the bulk. The long-time limit of the fluctuation 
function then follows from that of the retarded function using 

^ j d'^-'x ^'\{fiu,x),f{0)}) = (1 + 2nB(a;)) Im^'^Gl^u) (28) 

in frequency space. There is a similar relation for t^y correlators. 

3.7 Emergence of hydrodynamics 

The relevant Feynman diagrams are depicted in Fig. [31 Due to causality, the vertices 
will be restricted to lie inside the causal diamond emanating from the insertion points 
at the boundary, as depicted in Figs. [2] and HI 

In the r — t plane, the integration region Ai for the two vertices is the half-almond 
depicted in Fig. HI For the metric ([8]) its boundary can be described precisely. 



t = - ± (r^(r) - -), r*fr) = [tan ^(y)-tanh ^(y)] , y = [cosh(47rTr)] 2, 



T T ^ 



(29) 

where r is the (long) time separation between the two boundary insertions. For 
r ^ l/(27rT), the null geodesies are approximately horizontal lines, t ^ 0, t ^ t. 
In the r — 7- 00 limit, in which we are interested in, the tip of the half almond is 
approximated by 

ro(r) - ^e--^\ (30) 
Therefore, in this limit, the tip is exponentially close to the horizon. 
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Figure 4: Space-time structure of the integration region for the loop vertices ( "half- 
almond"). The external wavefunctions force the vertices to localize in the shaded 
region. Hydrodynamics interactions are generated in region A^> when both vertices 
are at radii of order T~^. 



The integration will not be dominated by the whole half-almond M. As will be 
seen in section (14.11) . the external wavefunctions decay exponentially on a microscopic 
time scale -C t as one moves vertically in the diagram. This will localize the 
integration regions to the shaded area in Fig. HI which has a thickness of order T"^. 
We will see this in more detail shortly. 

In addition, it will turn out that not all radii equally contribute. Heuristically, ac- 
cording to the scale-radius entry of the AdS / CFT dictionary, one can deduce that the 
radius at which a bulk vertex is evaluated corresponds to its non-locality scale in the 
boundary theory. In particular, vertices near the tip of the almond would correspond 
to interactions spread over a time interval ~ r, which may seem unexpected since 
vertices in the hydrodynamics are local in time and instantaneous. They correspond 
to a spread over a microscopic time interval ~ (see Fig. [6] below). With this in 
mind, we separate the bulk into two regions, Ai< and A^> at r = e, with e small. 
Basically e is chosen to separate the "tip" region of the half- almond from the rest. 
An adequate choice would be 



1 



~ -e--^^/\ (31) 



With this choice, the time separation between the vertices is large and order of r, 
everywhere in A^>. An important finding will be that the contribution from A^< is 
exponentially suppressed in r. The vertex integrations in region Aiy will give rise to 
the effective hydrodynamics interactions discussed in section |2l This will be shown in 
sections (14. 2 p and (14. 5p . This is how the one- loop boundary hydrodynamics emerges 
from the bulk gravity. 
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4 Pre-requisites for computing the diagrams 



4.1 External Wavefunctions 

Let us define the zero-momentum gauge field wavefunction = il)i^{r)e~^'^^ . The 
boundary to bulk propagator for and k'^x are just the wavefunction for these bulk 
fluctuations satisfying tp^ir) = 1 when r — )■ cxd. The appropriate boundary conditions 
are retarded for the field sourced by the field insertion at r = and advanced for 
the field sourced by the other one at time t = r@ The retarded boundary condition 
implies ipuji^) oc r~^2^ as r — 0. 

The linearized equations of motion for A^ give 

],dr[Gdri^^{r)] + \g''\uj^i;^{r) = 0, (32) 
y 

where Q = \/—gg^^. The graviton field h^^ can be discussed similarly instead Q = 
^—g. The wavefunction ■?/'(r) with these boundary conditions is 

^ ^^^^^ ^^^^^^ 

where tn = and 

271" J 

, (i + i) (-1 + ^ r(i-itt)) , , 

The general features of the time domain wavefunctio n fl33|) are easy to obtain. 
Note that the wavefunction (!33|) has quasinormal poles at lJ 

tt)„ = n(±l -z), n = l,2,.... (35) 

This indicates that A^ fluctuations relax exponentially on a microscopic time scales 
as the quasinormal modes are gapped and for all n ImrD„ ^ 0. By construction, near 
the boundary the wavefunction is approximated by A^if) ~ 5(t). Near the horizon 
the dominant r-dependence is due to the r~*" factor coming from expanding tanh. 
From this we conclude that 

Af(r,t)^Af(t-r,(r)), (36) 

which indicates that a given field profile falls into the black hole following null 
geodesies. 

An important property of ip{r,t) is that its time integral J dt%l){r^t) = 1 is inde- 
pendent of r. This follows from the w — )■ limit of f l5^ : The only solution to fl52]) 
with a; = which is nonsingular at the horizon, is a constant. The same is true for 
h^^. These general properties are expected to hold for any AdS black hole. 



^^This can be derived within the Schwinger-Keldysh formahsni. See subsection 15. II 

^^Similar computations have been performed in numerous previous works. See for example [27] . 



14 



4.2 Radial and time integrations for vertices in A^> 

In this subsection, we work out the bulk gravity interaction vertices corresponding 
to the non-hnear hydrodynamics. We will show that the dominant vertices in the 
hydrodynamic limit can be written in such way that they only involve the bulk op- 
erators t^^{r) and j^(r) defined in the previous section. This makes the holographic 
map between the hydrodynamic diagrams contributing to the long-time tail and the 
bulk quantum gravity amplitudes most transparent. For the moment we only focus 
on the region A^>, where several approximations are applicable. 

4.2.1 Dominant fields in hydrodynamic regime 

First we must identify the bulk fiuctuations which are long-lived. We consider fields 
in Fourier space with four-momentum components {u, 0, 0, k) which are taken to be 
all small, compared to T. We look for narrow peaks in correlation functions. The 
radial fiow equation f l26p can be trusted in region A^> in this kinematic regime. 

For the gauge fields, the transverse components do not exhibit diffusive poles; they 
relax on microscopic scales so they can be ignored in the long-time analysis. The At 
and Az components mix but because of Gauss' law, 

wF"* = kF'-', (37) 

we have that At/A^ ~ k/u 2> 1 since u ~ fc^/T for a diffusive mode. Thus, At 
dominates long-time correlations!^ Stress tensor correlators in the shear channel 
behave similarly, with htx dominating over hzx- 

Sound excitations involve the fluctuations htt, htz, hzz and h^x + . . . + hzz (spatial 
trace) [2H]- They are all of the same order, since u/k ^ Cs = 0{1) for sound waves. 
Note that all diagonal components of the spatial metric hij are involved. In fact, the 
on-shell fluctuations of hzz are equal to those of hxx = hyy = ... up to subleading 
viscous corrections (which are further suppressed by more powers of momentum). 
The spin-2 fluctuations can be neglected. This will be important later on. It can be 
shown as follows. The radial flow equation fl2^ for the spin-2 part of t^^ is initialized 
by the boundary hydrodynamic data. The initial data is provided by the linearized 
constitutive relation for the stress tensor 5tbdy — ^P^^'' + V^'^^u^^ + C(V^). The 
observation is that the spin-2 part is a viscous effect and therefore, suppressed by 
the momentum. Linearizing the deflnition ( 123|) . this demonstrates that the spin-2 
component of drh^j oc 5tbuik*j is also small in the bulk. One can evolve this data 
into the bulk using Since the boundary value of K'j is itself zero, its spin-2 

component will remain subdominant in the bulk. These estimates are outlined in 
Table [D 

4.2.2 Bulk interaction vertices: current 

Our flrst step is to couple the boundary operator j^, or equivalently the bulk pertur- 
bation 5 Ax sourced by it, to bulk three-point vertices. We concentrate on the "future" 

^°Note that this in only true in A^> away from the horizon. See subsection 14.51 
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Mode 


Frequency scale 


Fields turned on 


Charge diffusion poles 
Shear poles 
Sound poles 


u oc fcVT 
u oc P/T 
UJ (X k 


A > A, ~ {qlT)At 

htx > hzx ~ {q/T)htx 

htt ~ htz ~ hzz ~ /ixx = hyy (no spin-2) 



Table 1: Dominant bulk fluctuations in the long-term limit, with z being parallel to 
the spatial momentum. 



vertex near t = t. We begin with the gauge-gauge-graviton vertex. This vertex is a 
third variation of Yang-Mills' (or Maxwell's) action 

■T^pxM 



Va- 



d 



X 



^q^a^aSym 
5A^ 



d'^+^xSA^dMiSgdAi- 



-)] 



dAiiV^iSF^'^^h^ + SF^'^'h 



Notice the particular placement of the partial derivatives; this is a natural choice 
as it exhibits the boundary current as an integral over a "bulk current". When we 
estimate region A^< below, it is imperative to consistently keep up with the same 
convention, o 

In the region A^>, the individual terms appearing in the sum over indices in the 
second term are either /c-suppressed compared to other terms or vanish by our choice 
of gauge, so we safely ignore them. The loop fields obey Maxwell's equations, up to 
contact terms. Such terms would have the effect of pinching two vertices together, 
which is exponentially suppressed in A^> due to the external wavef unctions. It will 
not affect the long-time tail physics. In essence, a contact terms would turn two cubic 
vertices into a four-point vertex, which is negligible in A^> for the same reason. Thus 
the internal lines will be taken to be on-shell. Keeping only the contribution from 
M = r (the other indices being fc-suppressed), the remaining two terms reduce to 



A-A-g 



-9 



9l+i 



(39) 



According to Table [H fluctuations in Ftr dominate over those of Fir ~ {k/T)Ftr, 
which makes the second term in the bracket sub-leading. In any case, this term could 
not produce any long-time tail since it mixes a sound mode with a diffusion mode, 
which do not overlap in spacetime. Thus we keep only the first term with fi = t 

^5F*'' 



V, 



A-A-g 



d'^^^xbA, 



9l+i 



-drh'^t 



Linearizing the definitions (123bp and fl23ap this may be written as 



VA-A-g = levrG j d''+\6A,-£^ 



^""^Any other choice would be related to this one by a surface term. Notice that, by virtue of 
causality and smoothness of the external wavefunctions, and the boundary conditions at infinity, 
such surface terms associated to the integration-by-parts vanish. 
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The profile for SA^ has a time extent of order which is much smaller than the 
long time r, which separates the vertices. Since its time integral is 1, it can be 
approximated as a ^-function at the cost of 0{1/t) corrections 



VA-A-g 



d''-'x6f{T) 6t\{T) 
5j*(r)5t*.(r) 



drlQirGg^ 



+ 0(l/r) 



d 



e+p 



(40) 



the error being a relative error. The last line follows from the integral representation 
of the inverse momentum susceptibility 1 / [e+p) [22] . 

For theories with non-Abelian global symmetries the associated bulk gauge fields 
are non-Abelian, and we need to consider the associated gauge-gauge-gauge vertices. 
Keeping only terms with one radial derivative, using rotational invariance, one finds 
only products SA^drSA^ or 6A'^drSA^, which are sub-leading compared to PUI) in 
the hydrodynamic expansion. Therefore, these interactions are sub-leading (at zero 
chemical potential). 

In sum, in the region A^> we find that at zero spatial momentum 



6j\t) 5t*^(r) 
e+p 



+ 0{1/t) 



(41) 



A similar result holds for the "past" vertex near r = 0. 



4.2.3 Bulk interaction vertices: stress tensor t^y 

We repeat the same procedure for the bulk field h^^ dual to the boundary t^^ fluc- 
tuations. The corresponding bulk interaction is a graviton-graviton-graviton vertex. 
This vertex is a third variation of the Einstein-Hilbert action 



a-9-9 
1 

8^ 



6g6g6gSEH 



1 



d^+^xh^y5^- r^nn^i^ny^ 



167rG 



g{R + d{d-l)/e' 
I) ' 



2& 



(42) 



Similar to the case, the position of the derivatives in f l42|) is important. The double 
variation represents the insertion of the two gravitons running in the loop. Notice 
there is a factor of two coming from the fact that both h^y and hyx are sourced by 
the same profile. 

We evaluate the double variation in ( H2|) using Leibniz's rule. The quantity in 
the bracket is the equation of motion of the background and it vanishes when no 
variation acts upon it. When only one variation acts upon it this yields the equation 
of motion of the linearized perturbation, which vanishes up to contact terms. These 
are negligible in A^> as explained in the previous subsection. The equation (H2l) 
reduces to 



1 



9-9~9 



8nG 



Xy 



Rxy 



d{d-l) 
2£2 



'9xy 
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The last term manifestly has no second variation. Using equations of motion of the 
background plus those of the linearized perturbations, one can show that 6'^{g^^R) = 
g^uPR. Therefore, the second variation of the second term vanishes since there are 
no background off-diagonal metric components. We will have 



V, 



9-9-9 



8nG 



d 



d+l 



Xy 



-gh^y5'R,y. 



At this stage we can begin to make hydrodynamic approximations. This means we 
only keep terms in R^y which contain two radial derivatives 



g-9-g 



d 



d+l 



Xy 



-gh^'H'^ 



2K K"^ 



KK^y 



(43) 



The last term manifestly has no second variation. Following Table [H we keep only 
perturbations of diagonal components of the metric and of /i**. At this point the 
fact that long-lived spin-2 excitations are absent proves important; by rotational 
invariance, it is impossible to couple h^y to two diagonal modes. It is also impossible 
to have only one insertion because of parity in the index t. Therefore, only the 
terms whose index structure allow for two insertions are kept 



V, 



9-9-9 



^ j d''+^x^gh^y\gu\5K\5K\. 



Finally, utilizing J dt h^x ~ 1 plus the radial flow equations ( 126|) as above, we obtain 



9-9-9 



X St X St y 



dr , + C(l/r) 



-9\9 



e+p 



d'^~^xSt\ St' 



yi 



(44) 



the radial integral being as above. 

The graviton-gauge-gauge vertex is sub-leading. The reason is, by rotational in- 
variance (and keeping only terms with two radial derivatives), this can only couple 
to F^^F^"^ ~ j^jy which according to Tabled] is subdominant compared to (jH]). This 
implies that the graviton-gauge-gauge vertices do not contribute to long-time tails at 
zero chemical potential. From this, at zero spatial momentum 



^xy 



St'"" st'y 

e+p 



+ 0{1/t). 



(45) 



4.3 Bulk-Bulk propagators 

In the preceding subsection, we have ignored u and k dependent corrections in (1261) . 
This amounts to equating the bulk-bulk propagators of the currents in the loop to 
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their boundary limits. These boundary correlators can be found in [151 128] 



.2; 



U + I'^rfk^ 



fc^ {co + i^'jsk'^y — c^k"^' 

where for AT = 4 SYM theory, D= ^,E = p = |e = tt'^N^T^/S and r, r' are 

two arbitrary radii in Aiy. Note that the imaginary part of these retarded functions 
times 2T/uj reproduces the Fourier transform of the hydrodynamic expressions for 
the fluctuation functions ([3]). 

To consolidate this intuition and validity of the approximation ( H6l) . here we com- 
pute the bulk-bulk propagators in the hydrodynamic limit. As we will see in this 
section, even at small u and k expressions ( H6|) will break down at exponentially 
small r,r' in due to oscillations. Furthermore, if one applies the radial flow 
equations to a two-point function, the right-hand side of ( l26l) will generically con- 
tain a contact term. This leads to the appearance of contact terms oc 6{r — r') in 
f H6|) . This is generically expected to be the case for correlator with "momentum" 
variables. In sum, equation (B6|) holds at small u,k (or long times), where r and r' 
are not-exponentially-small. Also contact terms are expected in the expressions fH6|) . 

We only compute the retarded Green's function G^{r,r';p) for the operators j* 
explicitly. The Green's function is equal to the retarded field j*(r) induced by a small 
perturbation to the Hamiltonian 6H = — j*(r')e*^'^. That is to say G^{r,r';p) solves 
Maxwell's equation (with Ar = and = (— w, 0, 0, k) for definiteness) 

dr[^ ]+^k^ = J*, (47a) 



9d+i 9d+i 



•,2 



ikg''drA, + iu\g''\drAt = ^j=J"^ (47b) 



where the bulk current is 



.g\gtt\ 



9l+i 



9d+i 



with J* = 0. Note that the bulk current is conserved. By solving the second equation 
(]47b[) for drAz and substituting the result into a radial derivative of the first equation 
times C, we obtain a closed equation for j^{r;p) (and thus for (^^(r, r';p)) 
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where 



C 



9d+i 



The second derivative of a delta function can be removed by adding a contact 
term to the correlation functioro 



[dr{Cdr) + {\9''\^^-9''k')C] 



Gi{Ty-p) + 



6{r—r') 



-P(5(r-rO. (48) 



The retarded function G*^ obeys infalling boundary conditions at the horizon, 



Qtt ^ ^-iu^/2nT g^^^ QU ^ 



-iu)/2ttT jYoj g^j^^ Neumann condition (9j.G^ = at Vc- 



As seen from fl47ap . the Neumann condition follows from the Dirichlet condition on 
At and A^. 

The equation fj48|) is valid for any u, k. We now solve it in a hydrodynamical 
expansion. We first solve the differential equation (j^8|) with vanishing right-hand side 
following closely [15] and then patch the local solutions together to obtain Gr. The 
near horizon behaviour can be factored out by writing 

F(r) = e'-/"'^'-^(l + e), 

where e goes to zero at the horizon. This is the only regular solution. In the hydro- 
dynamic regime, e and its first derivative are uniformly small. The function e solves 
the equation 



(49) 



driCdre) = iudr (Cv^j + Cg''k^ + 0{ue, k^e). 
This can be integrated to give fl^e 

dre = iuj{^\ - i-) + 1 1 dr{Cg''k') + 0{ue, k'e) 



(50) 



where a is an integration constant. The regularity condition at the horizon implies 



a = lim,^o 1/(C^7^) B One obtains 



drF{r) 



Ca 



—iuo + 



0{uj'^ log -,uk\k^), 



(51) 



where S~^(r) = dr{Cg^^) has the interpretation of a "susceptibility" at scale r 
[22] . The wlog^ relative error signals the breakdown of ( ISTl) at exponentially small 
r due to oscillations. When u ^ 0, a second linearly independent solution is given 
by F*{r). From (!5T|) the combination which obeys Neumann boundary conditions at 

G{r) = (1 - z^)F(r) + (1 + t^)F*{r) 



Such contact terms are generically present in correlation functions of F^^ . In flat space, for 
instance, = — grf^^jTzj^ ~ —S{r—r')— g2^^2_fcii , which is structurally similar to P5|) . 
^■^This is the DC conductivity of the plasma. 
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where S = lim^^^oo S(rc). By gluing G{r) at large r to F{r) at small r and imposing 
), we finally obtain the bulk-to-bulk retarded function in the hydrodynamic limit 

, 6{r-r') _ _^^^G{r)F{r')e{r-r') + F{r)G{r')e{r'-r) 



C{r)g'' 2(u; + i|A;2 



Here the key observation is that in the hydrodynamic limit, the bulk-to-bulk prop- 
agator G^(r, r';p) loses its radial dependence (after removal of a trivial 5-function) 
and becomes identical to the boundary-to-boundary propagator. This can be traced 
back to the fc^ factor on the right-hand side of fHS]) and smallness of the right-hand 
side in (|26|). 

Equation (152|) is in perfect agreement with the expectation (l46l) with identifying 
D = (j/S (known as Einstein relation). We expect similar agreement for correlators 
of other fields. Thus, as promised at the beginning of this subsection, equations (l46ll 
hold up to contact terms and higher order uj and k corrections in the region A^>, 
where the radii are not exponentially small. Bulk-bulk propagators other than the 
retarded one, follow from the bulk KMS conditions f|T8l) . 



4.4 Gauge-fixing, Faddeev-Popov ghosts and causality 

At this point the reader might be puzzled that, while we have emphasized on a causal 
formulation of the problem, the gauge (fTOj) we have chosen is acausal. Causal gauges, 
such as light-cone, temporal or harmonic would have seemed more natural. On the 
other hand, the (spacelike) radial axial gauge is convenient so let us briefly pause to 
explain why this will not pose a real problem. 

In order to gauge-fix, we introduce Faddeev-Popov ghosts. Notice that the radial 
component of the metric fluctuations, hrr and h^^ are not fixed by the Symanzik 
procedure; this only fixes the induced metric on the boundary at r = Tc, while one 
needs to impose f lTOj) and f lTT]) everywhere including in the bulk. The gauge fixing 
action we take is 

Sg = j d'^+^x^ [X^^^n^huN + A n-A) , (53) 

where is the outward normal vector to the boundary. The fields A and A^^ are La- 
grange multipliers. The ghost action associated to this gauge fixing is easily obtained 

Sfp = j d^^^x^g [c*^n^ 5brst(^?m7v) + c 5^^si^{n-A)\ , (54) 

where 

5bkst{9mn) = ^mCn + VtvCm, 5brst(^-^) =n-d c + Am n-d c^, (55) 

where V is the covariant derivative with respect to the full metric. Note that the 
resulting ghost action is interacting and couples the gauge and gravitational ghosts. 
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The ghost inverse propagator in the gravity sector is simply given by 



N\-l 



6^S 



FP 



(56) 



Given that the ortho normal vector n has only radial component, it becomes clear that 
the ghost fields do not propagate in the field theory directions. A similar behavior 
is found in the gauge sector. At vanishing chemical potential the gauge sector is 
decoupled. 

A subtle issue here is that in order to invert (156|) , one needs to identify the bound- 
ary conditions for the ghost fieldsQ 

The unsettling issue seems to be that the ghosts itself can propagate on spacehke 
distances outside the causal diamond in Fig. [2J This would seemingly place our claim 
that only quadrant R is relevant to our formulation of the problem in danger. It is a 
well-known fact that the spacelike axial gauge suffers from such problems [30]. The 
gauge field propagators will also contain "longitudinal parts" that propagate along 
the radial direction and similarly violate causality. From standard arguments, when 
computing an on-shell amplitude like what we do here, the longitudinal modes will 
couple only among themselves (e.g., all propagators in a loop will be simultaneously 
longitudinal). We expect that this will produce a contribution of the same form which 
exactly cancels the acausal one from the ghosts in both regions A^< and A^>. Let 
us recall that the amplitude we compute is gauge invariant. One way to "regularize" 
the (spacelike) axial gauge (i.e., to make it causal) may be to obtain it only as a limit 
of a causal gauge, similar to the proposal of [29] . 

4.5 Contribution from region 

Near the horizon the metric takes the Rindler form 



The total volume of the half- almond (in the r — t plane) region A^< is of order e^. 
With the choice ( 13T|1 . the volume is exponentially small in r 



If the approximations which led to trivialized the radial integrations in ( HOl) and (jBl) 
did hold in the small-r region as well, computation of the vertices would have been 
already completed. This is so because these radial integrations would converge at 
their lower boundaries since the contribution from r < e scales like the volume (158|) . 

Some of the boundary conditions can be deduced from the BRST invariance of the resuhing 
gauge fixed action. For instance the BRST invariance implies that c'' = on the boundary. This 
can be intuitively understood; the presence of the regularized boundary at some r = breaks gauge 
invariance in the radial direction; having a non-zero on the boundary will change the boundary 
so we must set it to zero. In general, gauge-independence of physical observables should determine 
the boundary conditions. 



ds'^ = dr^ - {2TirTfdt^ + 5ijdx'dx^. 



(57) 
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However, we must be more cautious. For instance, we have neglected drAi compared 
to drAf, yet we know from Gauss' law that 

drA, = -l^^drA, ~ \^drA, (59) 

where z is taken parallel to the loop momentum (which may, or may not, be 
perpendicular to the external indices x and y). As seen from flS^ . the factor cu/fc^ gives 
a suppression away from the horizon which justifies neglecting current fluctuations 
over charge fluctuation. But in the near horizon region the suppression does not exist 
anymore due to diverging metric factors. The second field in the loop, i.e., g^^drh^ z 
in ( l39l) . satisfies an equation similar to (1591) . Using this, one can estimate the relative 
strength of a vertex involving current and stress in (l38ll (a term which was neglected 
in A^>) to that of involving charge and momentum fiuctuations 

d''+'x^SA,^drh\ ~ (-)2 / d^+^x—6A,jH\. (60) 

9 d-\- 1 z J \l 9 

This is logarithmically divergent at small r, since j* and t*^ remain order 1 at small r. 
The volume suppression flSSl) appears to be lost. In the next subsection we argue that 
despite this fact, the contribution stemming from region A1< is vanishingly small. 



4.5.1 One vertex in and one vertex in A^> 

In order to verify that the contribution from region A^< is small, we find it convenient 
to divide studying the radial integration for the one-loop amplitude, into two cases. 
First we concentrate on the case, where one of the vertices is located in A^< and 
the second resides in A^>. Suppose that the lower vertex at ri, lies in A1<, and the 
upper one at r2 is located in A^>. In A^< all fields propagate along null geodesies. 
Furthermore both fields in the loop must follow outgoing null rays in order to connect 
with the exterior vertex at r2. Since both momenta are approximately null and 
coUinear, one might expect that their Lorentz-invariant dot product is much smaller 
than their individual components. To be more precise, in the near horizon region 
using light-cone coordinates x^ = t ± r^,(r), one has 



g'" dM< 



>1<97V02 ^ g^ 9+0i9_02 + (1 ^ 2) < g''''dr^idr^2- 



(61) 



Let us now check whether the vertices indeed take this form. In order to do so, we 
need some estimates. In frequency space in the near- horizon limit, the possible radial 
profiles for the various fields can be readily found from the linearized equations of 
motion (for the metric perturbations, see [2H])- We find 
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Non-oscillatory modes in various channels are pure gauge modes. These pure gauge 
modes are generated by the residual gauge symmetry after gauge-fixings ( TTU]) and (ITT!) 
[28j . However, in section f l4.4p we have argued that the contribution from these pure 
gauge modes is canceled by ghosts loops, guaranteeing gauge invariance. Therefore, 
we should only worry about the "physical" part of the loop fiuctuations 

where n and u are any space-time indices. We start from expression (1381) for the 
gauge-gauge-graviton vertex 

VA-A-g = [ d'+'x ^d^[^{6F'^''h; + 5F^^hl - Uf^'^K)]. (63) 
J 9d+i ^ 

By the estimates (!62|) the only dangerous terms near the horizon, involve either two 
radial derivatives or two time derivatives contracted together with an inverse metric 
gtt _ '\Yg (^.g^^ neglect any explicit appearance of gu so we drop any field with a lower t 
index. Such terms and terms containing other type of derivatives are represented by 
(9(r°) inside the bracket. Using this it is easy to see that middle term in the bracket is 
sub-leading. Going on-shell as before, in the near horizon limit, the other two terms 
read 



-/^^(?^^(9mM,)(9^/iJ) - ^/^^(9MM.)(9^/ii) + 0(r°) 



VA-A^g^U'^'x^4^^g 
J 9d+i 

(64) 

With the same power-counting we find that the three-graviton vertex near the horizon 
is 

V,.,., = --^ f d'+'xV^h^y5'R,y (65) 



SttG 



8^ / ^^'^V^f' 



'/"'(OMh\)(d,ih„) - </"'\(dMhry)(BMK,) + 0{r'>) 



We see that both vertices (IMIl and (IMil enjoy the near- horizon cancelation 
they vanish in the eikonal limit. In order to estimate the sub- leading terms in the 
eikonal limit the following observation is important. We can estimate the sub- leading 
corrections to the null-geodesic approximation of fl62|) . For a typical loop fiuctuation 
(pi in the near horizon limit, we have 0j — )■ r^*'^(l -|- rii{r,oj, k)), where rj ~ r"^' with 
Re7j > in order to respect ( l62i) . The correction rj is computable perturbatively. If 
A = min{Re7j}, from equations ( 16^ and ( l66l) . one finds that the vertices are at least 
as small as 

VA-A-g,V,^g-g - [eit)]^''-"-'. (66) 

The integration over them in A^< is thus suppressed by e'^ which is exponentially 
small in TO 



Since the integrand of a vertex is not gauge- invariant by itself, for this argument to hold it is 
very important that the derivatives have been organized not to act on the external SAx or h^^ fields, 



exactly as was mentioned in section (|4.2 
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4.5.2 Both vertices in Ai 



There remains the case where both vertices reside in A^<, where r < e(r). The region 
M < covers the "tip neighborhood" of the half-almond. We divide the discussion 
according to the spatial momentum circulating in the loop, small k < T ot large 

Consider a fluctuation with small spatial momentum k running in the loop. Since 
this region is fiat, the relevant fields will follow either incoming or outgoing null 
geodesies. Most importantly, there are no correlations between the two types. This 
is because correlations can only arise when the wave "bounces" against the geometry 
at large r which takes a very long time (of order r). Pick an incoming loop fluctua- 
tion. When it scatters against the incoming wave {SA^ or hxy), it will remain purely 
incoming and un-deflected. Now when it is inserted into the second vertex, the cance- 
lation mechanism discussed previously will be effective and gives rise to the advertised 
suppression. Recall that since both vertices are constrained to the region A^<, sup- 
pression of only one vertex is sufficient. Diagrams with four-point vertices also exist 
in this region. They enjoy the same suppression. This is because such a vertex will 
connect either to three ingoing fields and one outgoing field or to three outgoing fields 
and one ingoing field. Since any contribution at large r is exponentially small in r 
(due to the external wavefunctions) , one can freely employ integration-by-parts to 
make all derivatives to act on the fields moving in the same direction^ 

There remains only the large k/T region. Manifestly this has little to do with 
hydrodynamics. Since these modes decay very fast in time at large r, their contri- 
bution is small asymptotically so we can freely perform integration-by-parts in them. 
Since at large k the exterior, curved region is separated by a large potential barrier, 
we may think of this contribution as a gauge-invariant one-loop scattering amplitude 
happening in just flat space. This is an amplitude for a purely infalling perturbation 
to reflect back from the horizon into an outgoing fluctuation. Note that, because the 
energies in the loop are exponentially large, this may not necessarily be computable 
within gravity; a stringy computation might be necessary. However, regardless of 
how it is computed, it appears extremely unlikely that such an amplitude could be 
anything but exponentially small (in r) - a black hole event horizon is not supposed 
to reflect perturbations even at quantum level. A simple check of this claim is that 
it is true for any local counterterm. In sum, although unlikely to be computable in 
pure gravity, it appears safe to ignore the contribution from region A^<. 

This completes our argument that with a natural choice for the action of the 
derivatives in section (14.21) . only the region A^> is important. 

large derivatives of the metric will not be generated in this process, because the Rindler metric 
is flat. 
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5 Assembling parts 



5.1 Schwinger-Keldysh Basis 

It is most convenient to make a detour by first switching basis, within the two- 
dimensional vector space in which the real-time fields 0i^2 "live", to the so called 
Keldysh (or (ra))0 basis. The Keldysh basis fields consist of an averaged (or retarded) 
field 0r and a differenced (or advanced) field 4>a 



where we have used the equilibrium relations (!T5|l for the symmetric function Grr = 



This basis is almost always the most efficient one to use. Taking into account the 
fact that Gaa = 0, the only contributing configuration of the Keldysh indices for a 
one- loop retarded amplitude is shown in Fig. El^a). All interactions vertices carry an 
odd number of (pa fields and we only need the vertices with exactly one a index. 

The physical meaning of this diagram should be quite transparent. It amounts to 
expanding a retarded two-point function to quadratic order in a background fluctua- 
tion, which is then averaged over an ensemble with variance Grr (the doubly-dashed 
propagator in the figure). This gives the full quantum-mechanical amplitude at one- 
loop. The "causality" is manifest in from Fig. |5t^a); the products of retarded prop- 
agators occur in such a way that both vertices must lie inside the causal diamond 
emanating from the insertion points of the external operators. This is how the black 
hole interior region is manifestly avoided, as previously advertised. 

5.2 Long-time tails 

From the effective vertices fHTj) and fj45l) and propagators fH6l) . the agreement between 
the bulk computation and the boundary hydrodynamics prediction is already obvious, 
since the radial coordinate has now disappeared from the problem. The remaining 
steps of the computation are nonetheless tedious. We combine fllT]) . (jlH]) and the 
diagrammatic results of the preceding subsection to obtain the one-loop retarded 
correlator for j^. At this stage we can safely go back to frequency space, since the 
radial integrations have been taken care of, alleviating the need of manifest space-time 
causality: 



(Pr = + h), 0a = 01-02- 

In this basis, the propagator is the 2x2 matrix 

/ Grr Gra \ _ f -^(G'r-G'a) (|+r2B(w 





(67) 



(68) 



^{G>+G<). 




rim G 



{co', k) 



G'^''''{co-Lo',k) + G'^{co',k) 



ImGj^' ^(tu—Lo', k) 



(69) 



u—u' 



Here "r" meaning "retarded" not to be confused with the radial coordinate. 
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r 



c 



(C) 



Figure 5: (a) Keldysh index structure of the one- loop diagrams contributing to a 
retarded two-point function, (b) A diagram which vanishes due to a closed loop of 
retarded propagators, (c) A diagram with a four-point vertex (which is subleading 
in the long-time limit). Arrows indicate the flow of time along retarded propagators. 
The cut propagator is the G^r propagator ( "symmetric propagator" ) . 



where w = e + p. We have used Grr — (2T/cj)ImGR in the loop, which is valid at 
small frequencies uj,u)' '^T. The integrations are straightforward but slightly tedious 

111 / . " 



1 



{u'+iDk^)[{uj-u'y + i^nk^)^] {uj'^ + {Dk^f){u - uj' + i-i^k^) _ 
TED-f, [D+^r, kl+kl 



w Jk D-f.^ uj + i{D+-fr,)k'^ 



d-l 
2 



When a; < we take the complex conjugate. Note on the second line we have dropped 
the sound contribution from G^'*^, which does not contribute at large times since its 
pole does not pinch with the diffusion one in [8j. Similarly, on the last line we 
have kept only the term which is non- analytic at cu = 0; this originates from small k 
and is not affected by the fact that the k integration is ultraviolet-divergent. 

The approximations made in this paper make it impossible to discuss the low- 
frequency regime of G, only the long-time regime, e.g. the non-analytic behavior at 
small u. 

Taking the imaginary part of ( I70|) and multiplying it by 2T/u) gives the fluctuation 
function. We write it in time domain 



T^E d-2 „A-d^ . ,nl-d duj 
w (47r(D+7^))^ ^ 2 ' '22 V-oo 2vr* 



^^)GTAt) = ^ ,_, r(^)sin(^^) / 



T^Ed-2 1 1 X 

w c^- 1 [47r(D+7^^)|t|]— t— +1 

which exactly reproduces the hydrodynamic prediction (!6all . Note that this is the 
same as J d^~^x ^^\{j^{t,x),j^{0)}). 
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For t^^ the calculation is similar but it uses the effective vertex ( H5|) . t^^ ~ 
t^Hy^/{e+p). By now it is easy to commute the use of fluctuation-dissipation rela- 
tions with the Fourier transform^, making the equivalence with the hydrodynamics 
computation in section ([6]), even more clear 

(^)G'™(t,A;=0) = [ [G'tf"(t,A;)G'*f^(t,A;) + G'tf^(t,A;)G*f"(t,fc)] 

Jk 

= Eq. (I6b]). (72) 

The bulk symmetric functions Grr in ( l46ll appear here. They are precisely equal to 
the boundary ones (E])- 



6 Concluding remarks 

So far, in the context of AdS/CFT, the phenomena associated to the non-linear hy- 
drodynamics, had mostly remained unexplored. In this presentation, we have shown 
how to reproduce hydrodynamic long-time tails generated by non-linearities from a 
one-loop AdS gravity computation. Universality of the hydrodynamic description was 
clearly observed within gravity computation. Only very basic entries of the AdS / CFT 
dictionary were used. 

Based on scale/radius relation in AdS/CFT, one might have expected that the 
long-time tails originate from the region in the immediate vicinity of the black hole 
event horizon, where the gravitational red-shift factors are large. In our setup, it 
turns out that this expectation is not entirely true. It is true that the very existence 
of the long-lived diffusive fluctuations is intimately tied to the near horizon region, 
however we find that the scale at which these fluctuations interact, is set by another 
scale. Indeed, the interaction vertices, from a bulk point of view, localize at the same 
radius at which susceptibilities receive most of their contributions (see the integral 
(]40|)), i.e., the scale T. 

The way to reconcile this finding with the scale/radius correspondence is to observe 
that the interactions are associated with the microscopic scale T and not the long- 
time scale r. In fact, in hydrodynamics, the interactions are exactly instantaneous 
(see Fig. E]). If the interactions had been localized in the near horizon region, they 
would have corresponded to nonlocal (in time) interactions from the boundary theory 
viewpoint. This could not have been reconciled with the hydrodynamics predictions. 

From a subset of our computation, contained in subsection 14. 2[ one can deduce 
the interaction between three boundary hydrodynamic modes. The interactions found 
from the bulk gravity are in agreement with the hydrodynamic prediction. Related 
computations have appeared before, to all orders in nonlinearity in hydrodynamics [6j . 
Our method is complimentary. We have shown that the leading bulk gravity Feynman 
diagrams with three external legs reproduce the leading non-linear interactions in 
hydrodynamics. 

Equation ([72]) is only exact in the soft approximation ne ~ T/u. In general, with obvious 
notation, the correct formula for a (rr) self-energy is Hrr = / ^{G^G^ + G^G^), where G^'^ are 
Wightman functions. 
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(a) 



(b) 



Figure 6: (a) N=4 SYM microscopic processes contributing to (b) a hydro dynamic 
vertex. In the regime considered in this paper, (b) is also equivalent to a gravity 
vertex. 

As expected, the "correct" degrees of freedom to use in the bulk are the (radial) 
momenta conjugate to the metric and gauge field. This is nice because the boundary 
limit of the (radial) conjugate momenta are exactly the boundary charges and currents 
which are "the variables" in hydrodynamics. These variables describe fluctuations of 
the horizon. In the hydrodynamic regime, they obey simple radial flow equations 

The one-loop computation appearing in this paper is technically novel, given its 
dynamical and real-time nature, although a different one-loop AdS computation, in 
the context of AdS/condensed matter physics, has appeared prior to this work in 

It is not evident from our computation whether the asymptotically AdS region is 
necessary for the long-time tail phenomenon to occur. It would be interesting to look 
for similar phenomena close to the horizon of black holes in asymptotically flat space, 
in suitable correlators and on suitable time scales. 

In the presence of continuous phase transitions, hydrodynamic description must 
include the order parameter itself as a degree of freedom surviving the long- wavelength 
limit. Including its effect might provide an interesting extension of this work. 
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